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Annotation  .. 


In  this  work  we  examine  the  statement  and 
numerical  solution  of  the  probl.em  of  the  develop¬ 
ment  of  ionization  instability  in  a  low-tempera¬ 
ture  magnetized  plasma.  Basic  attention  is 
devoted  to  the  dynamics  of  the  development  of 
instability  and  a  qualitative  analysis  of  the 
phenomenon.  W& "STlOtT'^the  effectiveness  of  the  ^ 
physical  model. used  as  the  basis  of  the  examination- 
As  a  result  of  the  numerical  calculations. we-  \ 

.wrfine  the  association  between  the  average 
parameters  of  the  plasma  at  process  3tages  close 
to  turbulence-  pvr  >Jf 
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CHAPTER  1 

INTRODUCTION 

Recently  much  attention  has  been  paid  to  the  3tudy  of  the 
phenomenon  of  ionisation  instability  in  a  nonuniformly  ionized  low- 
‘temperature  plasma  In  a  magnetic  field. 

First,  it  appeared  that  this  phenomenon  is  of  practical 
significance  for  MHD  generators  and  accelerators  using  atomic  gases. 

It  appears  simultaneously  wish  the  separation  of  electron  temperature 
from  the  temperature  of  heavy  particles.  In  a  sufficiently  strong 
magnetic  field  the  ionization  Instability  determines  the  effective 
conductivity  of  the  plasma. 

Second,  ionization  instability  apparently  leads  to  one  of  the 
few  forms  of  turbulence  which  is  most  accessible  to  theoretical  and 
experimental  Investigation.  The  relatively  simple  structure  of  the 
energy  equation  makes  it  possible,  in  a  number  of  cases,  to  construct 
simplified  mathematical  models  of  the  phenomenon \  we  are  able  to 
examine  an  already  developed  process,  and  not  merely  its  initial 
stage.  The  results  of  the  theory  describe  qualitatively,  to  a  certain 
degree,  the  hasic  regularities  known  from  experiment.  In  particular, 
in  a  number  of  cases  the  Increase  In  effective  resistance  with 
magnetic  field  is  predicted  theoretically  and  confirmed  experimentally. 
From  dimensional  concepts  it  Is  to  be  expected  that  in  a  turbulent 
plasma  with  magnetized  electrons  the  resistance  should  be  defined  by 
the  relationship 
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where  B  is  the  magnetic  field,  c  is  the  speed  of  light,  e  is  the 
electron  charge,  and  ne  is  the  electron  concentration.  This  formula, 
as  can  be  easily  seen,  corresponds  to  the  familiar  Bohm  formula  for 
turbulent  diffusion  of  plasma. 


» 

An  essential  understanding  of  the  phenomenon  was  also  achieved' 
in  experimental  study.  This  was  aided  by  the  short  duration  of  the 
characteristic  time  of  development  of  the  instability,  which  made  it 
possible  to  isolate  it  in  "pure”  form,  i.e.,  to  avoid  the  influence 
of  other  instabilities  (ion  sound,  etc.).  This  is  an  additional 
reason  for  Interest  in  ionization  instability,  whose  study  is  thus 
important  both  from  an  applied  as  well  as  a  general  physical  point  of 
view.  Up  to  the  present  time  a  great  many  experimental  and  theoretical 
works  have  been  devoted  to  this  phenomenon.  Let  us  give  a  brief 
summary  of  certain  of  these. 


Attention  was  first  paid  to  the  possibility  of  the  origination 
of  ionization  instability  in  [1]  and  [2],  and  also  independently  in 

£33. 

<*  */ 

Reference  [2]  gave  a  one-dimensional  nonlinear  theory.  ,  V* 


It  was  shown  that  the  development  of  one-dimensional  fluctuations;, 
leads  to  a  state  with  the  effective  resistance 


p 


cB 

2en  ’ 
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(2) 


The  quasi-linear  method  used  in  this  work  was  not  valid  far  from  the 
stability  boundary.  The  results  of  the  work  wers.  refined  by 
Zampalioni  [173  and  Solbes  [16J;  however,  the  final  result  of  [2j  was 
correct.  In  T 33  Ohm's  law  is  propos-cd  and  compared  with  experiment, 
taking  into  account  fluctuations  in  electron  and  gas  temperature.  The 
possibility  is  indicated  of  "complete"  instability  of  a  weakly-ionized 
plasma.  The  work  is  limited  to  the  frameworks  of  the  linear  theory. 
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One  of  the  first  experiments  is  described  in  [*i].  It  is  shown 
that  the  effective  conductivity  across  the  magnetic  field  decreases 
as  'Vg;  the  assumption  is  made  that  the  drop  in  conductivity  is 
explained  by  ionization  instability. 

In  [53,  using  an  electron-optical  converter,  there  was  obtained 
an  instantaneous  photograph  of  the  distribution  of  electron  density 
(and  current)  in  a  plane  perpendicular  to  the  magnetic  field.  There 
was  detected  the  existence  of  well-developed  striations  at  an  angle 
of  <v30°-il0o  to  the  current.  Qualitatively,  the  result  coincided  with 
theory,  although  a  complete  explanation  of  the  structure  was  obtained 
only  upon  numerical  solution  of  the  nonlinear  boundary  problem  [23-253. 
Further,  a  detailed  study  was  made  of  the  condition  of  origination 
of  instability,  the  influence  of  the  transport  of  resonance  radiation, 
.the  influence  of  current  along  the  magnetic  field,  and  systems  of 
conductors  creating  "rigid  equipotentials"  and  stabilizing  the 
instability.  The  basic  conclusion  of  the  work  was  that  suppression 
of  the  ionization  instability  is  almost  as  unsuitable,  energywise, 
as  reconciliation  with  its  existence.  Calculation  was  made  of  a  model 
of  nonlinear  development  of  ionization  instability  with  time,  which 
led  to  relationship  (2)  for  effective  resistance  of  the  plasma,  and  a 
model  of  developed  turbulence  was  proposed. 

Experimentally,  [6]  Is  a  development  of  [5],  It  gives  the 
results  of  a  photographic  study  of  instability  in  a  discharge.  The 
structure  of  the  nonuniformities  and  the  Influence  of  the  degree  of 
ionization  are  studied. 

Experimental  works  [7-10]  Investigate  the  development  of 
ionization  instability  in  a  disk  channel.  The  use  of  a  disk  channel 
makes  it  possible  to  avoid  the  influence  of  electrodes,  while  various 
types  of  pre-Ionicers  make  it  possible  to  study  the  dependence  of  the 
development  of  ionlzat-ion  instability  on  the  initial  conditions.  It 
is  shown  that  independently  of  conditions  at  the  inlet  ionization 
instability  occurs  and  there  is  saturation  of  the  effective  Hall 
parameter.  In  experiments  [10]  with  thermal  ionization  in  a  frozen 
flow  with  complete  ionization  of  the  admixture,  ionization  Instability 
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does  not  develop.  In  experiments  [93  it  was  detected  that  with  total 
ionization  of  its  cesium  admixture  in  helium, saturation  of  the  Hall 
parameter  appears  with  relatively  low  concentration  fluctuations.  In 
this  series  of  experiments  the  nature  of  the  saturation  of  the  Hall 
relationship  with  total  ionization  remains  unexplained.  In  work  [11] 
there  were  conducted  measurements  of  the  effective  parameter  in  an 
equilibrium  plasma  with  a  CC>2  additive.  It  was  shown  that  saturation 
of  the  Hall  parameter  is  nonexistent  up  to  values  ft  ^  5.  Thus, 
phenomena  in  the  disk  channel  actually  are  determined,  apparently, 
by  Instability. 

The  flow  of  inert  gases  with  alkali  metal  additives  in  a 
Faraday-type  MHD  channel  was  studied  in  [12-13].  It  was  shown  that 
at  high  current  densities  ionization  instability  has  a  strong  Influence 
on  the  effective  conductivity  and  the  Hall  parameter. 

Work  [l1!]  investigated  the  development  of  ionization  instability 
in  a  discharge  across  a  flow  and  a  magnetic  field.  Study  of  the 
geometry  of  the  fluctuations  showed  that  with  sufficiently  large 
current  densities  ionization  Instability  actually  developed. 

Work  [15]  theoretically  investigated  the  development,  with  time, 
of  a  nonlinear  one-dimensional  wave.  Further  calculations  are  given 
in  [16]. 

Theoretical  survey  [17]  discusses  in  detail  questions  of  the 
ionization  equilibrium  and  energy  balance  In  a  plasma  with  detachment 
of  the  electron  temperature  and  stability  of  a  nonequilibrium  plasma. 

Work  [18]  calculates  the  effective  conductivity  of  plasma  with 
two-dimensional  nonuniformity.  It  Is  shown  that  even  weak  isotropy 
of  the  fluctuations  leads  to  a  law  of  linear  increase  in  effective 
resistance  with  the  magnetic  field  for  a  given  amplitude  and  spectrum 
of  the  fluctuation. 

The  development  of  ionization  instability  in  a  system  with  given 
boundary  conditions  causes  basic  difficulties  even  in  linear  statement 
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of  the  problem  because  of  the  absence  of  natural  modes.  The  problem 
of  the  evolution  of  an  initial  perturbation  with  ideal  segmentation 
of  the  electrodes  has  been  solved  in  [193.  It  was  shown  that  initial 
perturbations  increase  nonexponentially ,  although  the  criterion  of 
stability  is  practically  the  same  as  in  an  unbounded  region. 


In  [20-21],  following  works  [2]  and  [53 ,  there  is  calculated  the 
quasi-linear  development  of  ionization  instability  near  the  threshold. 
It  is  shown  that  asymptotically  the  following  laws  are  satisfied: 


°eff  '  ncr/a 


ncr  *  °efr 


Later  these  relationships  were  used  as  the  basis  of  one-dimensional 
calculation  of  MHD  generators.  The  first  work  of  this  type  was  [22]. 


Works  [23-253  solved  numerically  the  problem  of  the  development 

of  ionization  instability  in  a  strongly  magnetized  plasma.  It  was 

found  that  the  initial  perturbation  develops  in  striations  which  then 

are  destroyed  and  the  picture  becomes  close  to  turbulent.  Finally, 

work  [26]  is  analogous  to  [23-253,  but  here  there  Is  examined  only 

the  Initial  stage  of  striation  formation. 

» 

On  the  basis  of  this  brief  survey  we  can  conclude  that  analytical 
examination  of  ionization  instability  is  extremely  difficult. 

Nonlinear  analysis  has  been  performed  only  in  one-dimensional 
approximation.  At  the  same  time  the  influence  of  boundary  effects, 
as  experiments  showed,  can  in  certain  cases  have  decisive  influence 
even  on  the  qualitative  nature  of  the  process.  However,  the  twc- 
dimensional  problem  Is  difficult  for  analytical  examination  even  in 
•linear  approximation. 

In  such  cases,  an  effective  method  of  investigation  can  be 
calculation  experiment,  which  helps  to  overcome  the  mathematical 
difficulties  in  the  theory  and  fulfills  the  physical  experiment  which 
gives  a  limited  amount  of  information  due  to  measurement  difficulty. 

This  present  work  is  a  development  of  [23-253  and  is  devoted  to 
a  discussion  of  a  number  of  numerical  experiments  on  the  development 
of  ionization  instability  in  various  models  of  nonequilibrium  plasma. 
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CHAPTER  2 


PHYSICAL  ASSUMPTIONS  AND  EQUATIONS 

We  examine  a  quasi-neutral,  dense,  low-temperature  plasma.  We 
will  consider  the  plasma  to  be  atomic,  i.e.,  we  will  exclude  from  our 
examination  the  interaction  of  electrons  with  molecular  compounds. 

Then  the  characteristic  time  for  energy  exchange  between  the  electrons 
is  considerably  less. than  with  heavy  particles  because  of  the 
considerable  difference  in  masses  (ions  and  neutral  atoms).  Because 
of  this,  the  energy  connection  of  the  electrons  with  the  heavy 
component  is  weak,  and  the  electron  gas  of  free  electrons  and 
electrons  occupying  discrete  energy  levels  of  the  atoms  forms  an  almost 
closed  equilibrium  system. 

This  Indicates  that  there  occurs  Maxwellian  distribution  of  free 
electrons  and  Boltzmann  distribution  of  electrons  which  are  in  the 
excited  levels  of  atoms,  characterized  by  electron  temperature  Te. 
Collisions  with  heavy  particles,  as  has  already  been  noted,  only  weakly 
influence  the  deviation  from  thermodynamic  equilibrium  of  the 
electrons,  but  define  the  value  of  electron  temperature  Tg. 

The  basic  reason  for  disruption  of  equilibrium  in  a  cold  plasma 
Is  the  radiation  of  the  electrons  of  excited  at?ms  in  a  discrete 
spectrum.  Here  a  correct  theoretical  description  of  the  radiation 
transport  process  is  given  only  by  the  kinetic  approach.  T^e  fact  of 
the  matter  is  that,  together  with  photons  In  the  center  of  the  line, 
having  small  free  path  lengths  because  of  resonance  absorption,  an 
essential  role  is  played  by  photons  on  the  "wings  of  the  line,'*  whose 
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frae  path  length  is  commens urate  with  the  dimension  of  the  plasma 
itself.  Radiation  leads  to  three  effects:  deviation  from  Maxwellian 
distribution  of  the  free  electrons,  deviation  from  Boltzmann 
distribution  of  electrons  on  excited  levels  of  the  atoms,  and  cooling 
of  the  plasma  as  a  whole.  Estimates  given  in  [171  together  with  the 
results  of  experimental  works  show  that  in  a  sufficiently  dense 
plasma  (n£  >  ncr  ^  10  1/cm  ;  ng  is  the  concentration  of  free 
electrons)  at  temperatures  Tg  0.3  eV  we  can  disregard  the  deviation 
of  the  electron  gas  from  equilibrium  because  of  radiation.  But  in 
dense  atomic  plasma  as  well  the  contribution  of  radiation  to  the 
electron  energy  balance  is  completely  commensurate  with  other  transport 
processes. 

Let  us  discuss  further  a  simple  model  of  an  atomic  two-temperature 

plasma.  The  light  component  consists  of  an  equilibrium  electron  gas 

with  temperature  Te,  concentration  of  free  electrons  ne,  and  electrons 

in  the  excited  level  of  the  atoms.  The  heavy  component  consists  of 

ions  and  neutral  atoms  with  concentration  n,  +  na  and  temperature  T. 

To  be  specific,  we  have  in  mind  an  inert  gas  (e.g.,  Ar)  with  a  slight 

addition  of  a  lightly  ionised  admixture  (e.g.,  Cs).  In  this  case  the 

concentration  of  electrons  and  ions  n^  coincides  with  the  number  of 

ionized  atoms  af  the  admixture  per  unit  volume  and  does  not  exceed 

the  concentration  of  the  admixture  n_  (we  consider  only  single 

s 

ionization) . 


With  application  of  an  electrical  and  magnetic  field,  there 
occurs  in  the  plasma  an  electric  current  whose  density  J  is  associated 
with  the  strength  of  the  electrical  $  and  magnetic  H  fields  by  the 
generalized  Ohm’s  law  [31-32].  In  our  case  we  can  disregard  Inertia 
of  the  electrons  and  "slip”  of  the  ions  relative  to  the  neutral 
atoms.  Conservation  of  momentum  of  the  electrons  for  a  Lorentz  gas 
under  these  assumptions  gives  Ohm’s  law  in  the  form 


J  +  o  x  ft  *  cr£’ - — Vp  ,  5  =  -SS- 

d  °  mv  ^e*  in  c 

e 


(2.1) 


where  a  is  the  conductivity,  defined  by  the  expression 
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(2.2) 


m  cv 
e 

Pe  is  the  gas  pressure  of  the  free  electrons;  is  the  electrical 
field  in  a  system  of  coordinates  moving  with  the  mean-mass  flow  rate. 

Here  and  henceforth  we  will  use  the  standard  designations : 

e  is  the  electron  charge, 
me  is  the  electron  mass, 
c  is  the  speed  of  light, 

M  is  the  mass  of  an  atom, 

I  is  the  ionization  potential, 

is  the  collision  frequency  of  an  electron  with  ions, 
v„  is  the  collision  frequency  of  an  electron  with  atoms, 

a 

k  is  the  Boltzmann  constant. 


Now  let  us  write  the  energy  balance  of  an  electron  gas  per  unit 
volume  per  unit  time.  According  to  the  first  law  of  thermodynamics 
it  can  be  represented  as  follows: 


^Cne(I  +  |kTe)]  -  Q  -  Qs  -  pe  x 
X  div  v  -  ng(I  +  |kTe)  div  v  -  div  q. 


(2.3) 


Let  us  examine,  in  sequence,  all  components  of  the  right  side  of  the 
equation.  The  first  term  defines  the  magnitude  of  &he  contribution 
of  Joule  heating  per  unit  volume; 


Q  * 


(2.10 


The  second  term  characterizes  the  output  of  electron  energy  during 
collisions  with  heavy  particles: 


-  ■V)efk<Te  -  T£  K  r  (2-5) 

a  ~  "a 

where  6a  is  a  constant  which  is  defined  by  mechanisms  of  collision  of 
an  electron  with  an  a-type  heavy  component. 
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The  work  of  electron  pressure  is  described  by  the  third  term  in 
equatior  (2.3)*  The  other  components,  in  sequence,  show  the  change 
in  energy  of  electrons  due  to  transport  processes:  convective 
radiation  and  heat  flux  q.  The  connection  between  concentration  ng 
and  temperature  Te  of  the  electrons  is  given  by  the  equation  of 
ionization  kinetics.  Let  us  assume  that  there  occurs  only  thermal 
ionization  by  electron  impact: 

a+e*i+e+e  (2.6) 

while  the  inverse  process  is  recombination  with  triple  collision: 

i+e+e»a+e.  (2.7) 

The  equation  of  the  kinetics  for  such  a  reaction  has  the  form 


dn 

-Ot  *  a*ne*na  “  3nine‘  (2‘8> 


The  constants  of  ionization  a  and  recombination  g  rates  are  connected 
by  the  principle  of  detailed  equilibrium: 

f*K(Te).  (2.9) 


The  equilibrium  constant  K(Tg)  in  an  equilibrium  electron  gas  is 
defined  by  the  Saha  formula: 


K(Te) 


nenl 

n_ 


2g1f2«lkTe1 


(2.10) 


where  g ^  and  g&  are  the  static  weights  of  an  ion  and  a  neutral  atom, 
respectively. 


The  description  of  this  model  of  the  plasma  is  closed  by  the 
Maxwell  equations.  As  usual,  let  us  disregard  the  displacement 

1  3  k  ll  <rr  4. 

currents  —  as  compared  with  conductivity  currents  — —  J.  Let  us 
Col/  c 

4  i|4r 

also  assume  that  the  induced  magnetic  field  h  — JqLq  (jQ  is  the 

average  value  of  electric  current  density,  Lq  is  the  characteristic 
dimension  of  the  plasma)  is  considerably  less  than  the  applied  field 
S.  In  this  case  the  equations  of  electrodynamics  express  the  law  of 
conservation  of  electric  charge: 
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(2.11) 


div  jf  *  0 

and  potential  oT  the  electric  field 

curl  1  =  0.  (2.12) 

Let  us  assume  that  the  external  magnetic  field  in  which  the  plasma 

« 

is  located  has  only  one  component.  Then  the  direction  along  the 
magnetic  field  £  is  isolated,  since  along  it  the  mobility  of  the 
electrons  is  considerably  higher  than  in  any  other  direction 
perpendicular  to  5.  Therefore  the  problem  can  be  considered  to  be 
two-dimensional,  and  it  can  be  examined  in  the  plane  normal  to  the 
magnetic  field.  Without  restricting  generality  let  us  assume  that  the 
plasma  occupies  a  rectangular  volume  ABCD  (Fig.  1)  whose  walls  are  an 
ideal  dielectric  excluding  electrodes  ab  and  cd,  which  are  either 
ideally  sectioned  or  compact  and  ideally  conducting.  On  the  dielectric 
portion,  obviously,  j  *  0..  On  an  ideally  sectioned  electrode  we  are 
given  the  electric  current  density  j  ,  while  on  the  ideally  conducting 
electrode  the  condition  ET  *  0  should  be  satisfied  and  we  are  given 
the  total  current  J  flowing  between  the  electrodes. 


Fig.  1. 


This  statement,  borrowed  from  an  experiment  in  which  ideally 
sectioned  electrodes  completely  occupy  the  horizontal  portions  of  the 
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boundary,  make  it  possible  to  trace  the  development  of  ionization 
instability  caused  by  a  perturbation  introduced  into  a  stationary 
field  of  electric  currents  [23-25].  When  on  the  horizontal  boundary 
there  are  portions  of  an  ideal  dielectric,  the  problem  simulates  the 
development  of  instability  in  an  MHD  channel.  Finally,  the  presence 
of  a  compact  ideally  conducting  electrode  makes  it  possible  to 
consider  its  influence  on  the  turbulent  plasma. 


CHAPTER  3 


DIMENSIONAL  ANALYSIS 

Let  us  formulate  more  specifically  the  statement  of  the  problem 
which  will  be  examined  in  what  follows.  Let  us  assume  that  in  a 
quiet  plasma  there  is  no  radiant  heat  transport  or  heat  flow.  We  will 
also  consider  that  the  electron  temperature  considerably  exceeds  the 
gas  temperature  Te  »  T. ,  but  at  the  same  time  the  average  thermal 
energy  of  the  electron  is  low  in  comparison  with  the  ionization 
potential  I  »  ~kTe.  The  collision  frequency  of  an  electron  with  the 
neutral  component  is  defined  by  the  expression 


va  *  Cogng  +  aa(ns  - 


while  the  frequency  of  electron  collisions  with  Ions 


vi  *  ai<l vei >ni* 


(3.1) 


(3.2) 


Here  o^,  as,  and  o ^  are  the  collision  cross  sections  of  an  electron 
with  atoms  of  an  inert  gas,  the  admixture,  and  ions,  respectively; 
<|ve|>  is  the  average  value  of  the  thermal  velocity  of  an  electron; 
rig,  ng  -  n^,  and  n^  are,  respectively,  the  concentrations  of  atoms 
of  an  inert  gas,  admixture,  and  ions.  The  total  frequency  of 
collisions  undergone  by  an  electron 


vi  +  v 


(3.3) 


We  will  assume,  in  addition,  the  following  dependence  of  cross  sections 

o.,  a0,  and  o  and  the  average  thermal  velocity  < | v _  j  >  on  temperature 
i  j  6  e  ' 


y 


(3.4) 


°i  *  °lTeC»  °s  *  a2»  ag  =  a3 


|vel>  '  feTe] 


(3.5) 


where  a.^,  a2,  and  are  constants.  Then  for  the  total  number  of 
collisions  we  get  the  expression 


v  -  A'neT;3/2  +  A'T^/2(ns  -  ne>  >  A -^/2ng 

A'  -  c  f3K] 1/2  A.  .  „  f^l1/2  a'  -  o  (2!l]1/2 

1  al  (m  J  *  a2  2[mJ  *  A3  3[m  J 


(3.6) 


Under  these  additional  assumptions  the  initial  system  of  equations 
obtained  in  the  previous  section  assumes  the  form; 


J  +  J  *  8  -  o2>  -  ^vPe;  a  -  -a,  8  .± 

I-jnr-  *  4r  ~  'Sn  v-j—kT  ,  v  *.Aim  T“^2  + 
dt  a  e  M  e*  lee 

+  A2Te/2(nS  -  V  +  , 


'e _  .2  a  2Slf2TOe!tTel3/2  "w« 


(3.7) 


(3.8) 


~dt  *  anane  “  6nine»  B 


3  ga  l  .  h 


div  J  *  0 
curl  £  *  0. 


(3.9) 


(3.10) 

(3.11) 


Let  us  reduce  (3.7)*-(3.H)  to  dimensionless  form.  As  the 
characteristic  scales  of  the  problem  we  select  the  Initial  value  of 
electron  concentration  nQ  at  a  certain  point  in  space,  temperature 
Tq,  and  electrical  conductivity  aQ.  The  scales  for  measuring  the 
other  magnitudes  are  expressed  from  equations  (3*7)— (3-9) •  From 
expressions  (3.7)  for  electrical  conductivity  we  have  the  scale  of 
measurement  of  the  collision  frequency; 
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(3.12) 


0  ~  a0me* 


The  dimensionless  expression  of  equation  (3.6)  assumes  the  form 


v  '  AlneTe3/2  +  A2Te/2(nE  '  na>  +  k?Y\ 


(3.13) 


Here  ng  is  the  reciprocal  of  the  degree  of  ionization  of  the  admixture 
at  the  initial  moment  of  time.  The  values  of  the  dimensionless 
constants  A-^,  A2,  and  A^  are  defined  from  the  condition  v  *  1  at  the 
initial  moment  of  time,  i.e., 


A,  +  (1  -  ne^A0  +  n„A0  *  1 
1  s  2  g  3 


(3-lJ0 


and  the  equalities 


A2  a2>  A3  7 3 


(3.1V) 


From  (3. 11!)  and  (3.1V)  for  A-^,  A2,  and  A^  we  get  the  values 


A1  *  ai  +  (ns  -  l)o2  +  ng03*  A2  “  ai  +  («s  -  D°2  +  ng®3* 


3  ax  +  (ns  -  l)o2  +  nga3* 


(3.15) 


From  energy  equation  (3.8)  it  follows  that  the  characteristic  time  of 
the  problem 


fc0  =  *  K  “  6M 


(3.16) 


whence  we  obtain  the  scale  of  the  measurement  for  electric  current 
density 


m  • 


(3.17) 


From  Ohm’s  law  (3-7)  the  scales  for  electrical  field  Eg  and  the  mean- 
mass  rate  Vg  have  the  values 
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E0  *  *V°0»  V0  *  CH^ 


(3.18) 


4 


while  the  measurement  scale  for  the  magnet4 c  field 


mcv 


0 


(3.19) 


is  such  that  the  dimensionless  value  of  the  magnetic  field  is  equal 

to  the  Hall  parameter  at  the  initial  moment  of  time.  Let  us  estimate 

the  contribution  of  electron  pressure  p  *  kn  T  to  Ohm’s  law  as 

e  e  e 

compared  with  the  other  terms.  Let  us  examine  the  relationship 


-  CPeO 
^  mevQL0^0 


(3.20) 


Here  pcq  is  the  characteristic  value  of  electron  pressure.  Using 
-(3.12)  and  (3.18)  we  get 


Y 


k*e0 

efeOL'o‘ 


(3-21) 


For  MHD  converters  and  most  experiments  the  following  condition  is 
valid : 

» 

|kTe  «  eEL 


i.e.,  the  mean  thermal  energy  of  an  electrode  is  low  as  compared  with 
the  characteristic  difference  in  potentials.  From  this  it  follows 
that  the  electron  pressure  gradient  in  Ohm’s  law  can  be  disregarded. 
Let  us  turn  to  the  equation  of  ionization  kinetics  (3.9).  With  low 
degrees  of  ionization  we  see  from  this  equation  that  the  character¬ 
istic  time  of  a  change  in  concentration  of  electrons  is  expressed 
by  the  value 


T 


1 

ne°0 


(3.22) 


where  aQ  is  the  scale  of  measurement  of  the  ionization  rate,  which  is 
expressed  by  the  formula 


a  *  /  vf(v)dv 

ne 

IW  >  v0 


(3.23) 


IT 


«a  is  the  interaction  cross  section  of  an  electron  with  an  atom,  vQ 
is  the  minimum  velocity  of  the  electron  of  ionizing  atoms  of  the 
admixture,  and  f(v)  is  the  distribution  function  for  free  electrons 
with  respect  to  velocities.  Considering  that  Maxwellian  velocity 
distribution  is  valid  for  electrons,  lor  a  we  get  the  expression 


a  ~  a 

a 


1/2  -I/kT 
e  e 


(3.24) 


respectively  for 


(3.25) 


Comparing  the  characteristic  times  tQ  and  x,  we  find,  assuming  that 

v0  -  W 


*kTe0^0 

w 


8*VeOk 

— ffi  ■ 


(3.26) 


Henceforth  we  will  assume  that  the  plasma  has  parameters  for  which 
x/t0  «  1*  i.e,,  the  thermal  ionization  is  quasi-equilibrium  and, 
consequently,  instead  of  equation  (3*9)  we  can  use  the  principle  of 
detailed  equilibrium: 


¥e 

n 


2Sif2™eMel3/2  "^e 

TT" 6  ' 


(3.27) 


In  dimensionless  form  (3.2?)  assumes  the  form 


n 


s 


n 


e 


(3.28) 


where 

kQ  »  I/kT0  (3.29) 

and 
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(4.3) 


«  nii. 

3y  ‘ax* 


Thus,  in  the  general  case,  for  equation  (4.2)  on  a  section  of  the 
boundary  vre  are  given  the  first  boundary-value  problem 

¥lr  *  v(s)  (4.4) 

where  s  is  a  length  of  arc  along  boundary  r,  while  in  the  remaining 
section  (it  may  be  absent)  we  have  the  condition  (4*3),  i.e.,  there 
occurs  a  mixed  boundary- value  problem. 


To  solve  the  equation 


dt 


<2 

-  n  vT 
o  e  e 


v  '  AlneTe3/2  +  A2(ns  -  ne)Te/<:  + 
we  are  required  to  give  only  the  initial  condition 

n|t«0  *  n0(x»  y)‘ 


(4.5) 


(4.6) 


The  form  of  the  function  nQ(x,  y)  will  be  selected  below  when 
examining  specific  examples.  As  linear  analysis  shows,  problem  (4.2) 
(4.6)  is  correct,  but  unstable.  For  the  unbounded  linear  problem  the 
perturbations  increase  exponentially  with  time,  while  for  the  linear 
boundary- value  problem  [19]  the  increase  in  perturbation  occurs 
nonexponentlally .  Let  us  note  further  that  since  the  differential 
operator  of  equation  (4.2)  is  not  self-conjugate,  in  it  actual  eigen¬ 
functions  are  absent. 


Nonlinear  system  (4. 2)- (4. 6)  was  Integrated  numerically.  In  this 
regard  let  us  note  that  there  are  no  diffusion  terms  in  differential 
equation  (4.5).  However  during  numerical  solution  of  the  problem 
there  appears  an  averaging  scale  X  equal  to  the  step  of  the  difference 
network  h.  When  solving  the  problem,  taking  into  account  diffusion 
processes  (thermal  conductivity,  radiation,  etc.),  the  diffusion 
length  of  these  processes  is  the  characteristic  dimension  of  averaging. 
From  this  standpoint  the  step  of  the  difference  scheme  can  be  treated 
as  the  diffusion  length  of  certain  physical  properties.  Thus  the 
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n#  * 


3/2 


(3.30) 


_1_  2gjLr2rrmekTe0| 

n0  2sal  ^ 

Bearing  in  mind  that  the  scales  n0  and  TV  are  their  initial  values, 
satisfying  (3.27),  the  magnitude  can  be  found  from  the  equality 

k0 

n*  =  e  / (ng  -  1).  (3. 3D 

Finally,  let  us  formulate  the  equations  which  describe  the  process 
in  dimensionless  variables: 


dn 


— £  si - n  •  v*T 

dt  o  ne  v  le 


v  «  A^t;3'2  +  Aj(ns  -  %»l/2  +  A3ngT^/2 


(3-32) 

(3-33) 


n. 


ns  ‘  ne 


„,T3/2e‘k<>/Te 


J  J  x  $  *  ff(^  +  v  x  $),  &  «  S/v,  a  *  n/v 

div  J  »  0 
curl  2  «  0. 


(3.3*0 

(3.35) 

(3.36) 
(3- 37) 


Let  us  examine,  as  an  example,  values  of  uhe  dimensionless 
parameters  for  a  plasma  of  cesium  vapors.  Let  us  assume  that  at  the 
initial  moment  t  «  0  the  electron  temperature  T  -  2900°K  and  the 
degree  of  ionization  a  «  1/3  (the  electron  concentration  nfi  *  6.9*10 
1/cm^).  Then  we  find  that  kQ  *  10  and  ng  *  3.  Knowing  the  conductivity 
of  such  a  plasma  at  the  initial  moment  of  time  oQ  we  can  find  the 
values  of  the  remaining  dimensional  parameters. 
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CHAPTER  4 


MATHEMATICAL  STATEMENT  OP  THE  PROBLEM 


The  distribution  of  electric  fields  E„,  E„  and  currents  J  . 

x  y  x  y 

is  conveniently  described  using  the  vector-potential  of  the  electrical 
current  $.  In  the  examined  case  of  a  two-dimensional  problem  we  can 
consider  that  only  its  z-component  is  different  from  zero;  this  we 
designate  by  ip.  By  definition  of  the  vector-potential: 


.  i 4  i  ..2! 

*  3y*  Jy  3x* 


(4.1) 


Thus,  knowing  function  ip,  we  can  seek  the  distribution  of  currents 
Jx,  Jy-  Excluding  Jx,  jy  and  Ex,  Ey  from  equations  (3-35)~(3-37) , 
using  (4.1)  we  get  an  equation  for  ip: 


JLfi  it]  +  JJl  it]  +  JLfo  it]  .  JLfo  it]  ,  0 

3x(o  3xJ  3y [a  3yJ  r  3xja  3yJ  3y|a  3xJ  u 


(4.2) 


£  *  1/v,  a  *  n/v. 


As  has  already  been  mentioned,  we  will  examine  rectangular  region 
ABCD  (Pig.  1),  whose  boundary  is  an  ideal  dielectric,  excluding 
electrodes  ab  and  cu.  On  the  dielectric  sections  aACc  and  bBBd, 
ij>  =  const.  Considering  the  arbitrariness  in  defining  ip,  let  us  set 
ip  *  0  on  aACc,  and  ^  =  JQ  on  bBDd,  where  JQ  is  the  current 

passing  between  the  electrodes.  Electrodes  ab  and  cd  will  be 
considered  either  ideally  sectioned  with  constant  current  density 
J  =  Ja/Lab>  or  ideally  conducted.  In  the  first  case,  on  sections 
ab  and  cd  we  are  given  function  ip  as  linearly  Increasing  from  0  to  tj/^. 
In  the  second  case  E  =  Ex  =  0  on  sections  ab  and  cd  leads  to  the 
need  for  giving  a  condition  with  the  "directional  derivative" 
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Thus,  in  the  general  case,  for  equation  (4,2)  on  a  section  of  the 
boundary  we  are  given  the  first  boundary-value  problem 

^|r  =  v(s)  (4.4) 

where  s  is  a  length  of  arc  along  boundary  r,  while  in  the  remaining 
section  (it  may  be  absent)  we  have  the  condition  (4.3),  i.e.,  there 
occurs  a  mixed  boundary-value  problem. 


To  solve  the  equation 


n  vT 
e 


e 


WeV2  +  A2<"S  -  «e>Te/2  + 


A-n  T 
3  g  e 


1/2 


(4.5) 


we  are  required  to  give  only  the  initial  condition 

n|t«0  “  no(x>  y)*  C^‘6) 

The  form  of  the  function  nQ(x,  y)  will  be  selected  below  when 
examining  specific  examples.  As  linear  analysis  shows,  problem  (4.2)- 
(4.6)  is  correct,  but  unstable.  For  the  unbounded  linear  problem  the 
perturbations  increase  exponentially  with  tlm^,  while  for  the  linear 
boundary- value  problem  [193  the  increase  in  perturbation  occurs 
nonexponentially .  Let  us  note  further  that  since  the  differential 
operator  of  equation  (4.2)  is  not  self-conjugate,  in  it  actual  eigen¬ 
functions  are  absent. 


Nonlinear  system  (4. 2)- (4. 6)  was  integrated  numerically.  In  this 
regard  let  us  note  that  there  are  no  diffusion  terms  in  differential 
equation  (4.5).  However  during  numerical  solution  of  the  problem 
there  appears  an  averaging  scale  X  equal  to  the  step  of  the  difference 
network  h.  When  solving  the  problem,  taking  into  account  diffusion 
processes  (thermal  conductivity,  radiation,  etc.),  the  diffusion 
length  of  these  processes  is  the  characteristic  dimension  of  averaging. 
Prom  this  standpoint  the  step  of  the  difference  scheme  can  be  treated 
as  the  diffusion  length  of  certain  physical  properties.  Thus  the 


change  in  step  in  the  difference  net  leads  to  a  change  in  the 
characteristic  scale  of  the  problem.  Nonlinear  system  (4.2)-(4.6) 
was  numei'ically  integrated  as  follows.  First  we  found  the  solution 
of  the  steady-state  equation  (4.2)  with  coefficients  calculated  from 
the  values  of  electron  concentration  ne  from  the  previous  time  level. 
Then,  according  to  the  difference  analog,  equations^  (4.5)  were 
recalculated  for  values  of  ng  at  the  new  moment  of  time.  The  boundary 
value  problem  for  equation  (4.2)  has  a  number  of  specific  features 
which  should  be  taken  into  account  when  solving  it  numerically.  The 
basic  ones  of  these  are  as  follows: 

1.  The  values  n,  v,  and  a  can  be  highly  variable  in  space. 

2.  The  differential  operator  is  not  self-conjugate. 

3-  On  an  ideally  conducting  electrode  there  is  a  condition  with 
a  "directional  derivative." 

4.  The  value  a  on  a  certain  previously  unknown  section  can 
vanish. 

The  calculating  algorithm  was  developed  on  the  basis  of  [27-28}. 

A  detailed  exposition  of  the  methodological  questions  is  given  in  [29]. 


CHAPTER  5 


RESULTS  OF  NUMERICAL  CALCULATIONS.  DISCUSSION 

1.  Let  us  begin  the  examination  with  the  simplest  statement  in 

which,  nonetheless,  the  characteristic  features  of  the  studied 

phenomenon  are  expounded.  Let  us  assume  that  in  rectangular  region 

ABCD  (Fig.  1)  there  is  a  plasma  with  constant  concentration  of  free 

electrons  n  *  1. 
e 

The  external  magnetic  field  is  also  constant  throughout  the 
region.  In  the  plasma,  between  the  ideally  sectioned  electrodes 
ab  ■  AB  »  1  .and  cd  »  CD  *  1,  there  flows  a  constant  and  uniform 
electric  current: 

Jx  "  °»  Jy  =  _1- 

Walls  AD  =  2  and  BC  *  2  are  made  of  an  ideal  dielectric  on  which 
i  -  0.  The  assumed  values  of  the  electric  current  density  and 
electron  concentration  are,  formally,  the  solution  of  steady-state 
equations  (4. 1-4. 6)  which,  however,  as  will  be  seen  from  what  follows, 
is  unstable. 

Let  us  assume  that  at  moment  t  =  0  there  is  a  perturbation  of  the 

three-dimensional  distribution  of  electron  concentration.  We  will  be 

interested  in  the  further  time  behavior  of  concentration  n  ,  of  the 

vector  field  of  the  electric  current  density  j  ,  j  ,  and  also  the 

a  y 

Integral  and  average  characteristics.  These  latter  play  a  very 
substantial  role  in  the  study  of  the  process.  In  this  present  work 
we  Introduce  into  examination 
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<n> 

l 

// 

n.dxdy 

bABCD 

ABCD 

<a> 

*  1 

If 

rjdxdy 

SABCD 

ABCD 

<Ex> 

1 

// 

E„dxdy 

sabcd 

ABCD 

A 

<E  > 

IS 

E„dxdy . 

y 

sabcd 

ABCD 

y 

When  constructing  the  average  <E  >  and  <E  >  it  was  considered  that  the 

x  y 

region  in  which  the  electric  field  is  effectively  present  comprises 

part  of  region  ABCD.  To  be  specific,  we  will  arbitrarily  consider 

that  this  is  the  Interelectrode  area  sa^cd‘  a<^ifci°n  we  will 

trace  the  minimum  n  and  maximum  n  _  concentrations  of  free 

min  max 

electrons  in  region  ABCD.  Of  great  interest  are  the  effective 
conductivity  of  the  plasma 

°eff  *  <3>/<B3> 
and  the  effective  Hall  ratio 


In  the  examined  statement 


Beff 

<i>  1 

aeff 

geff 


-  <ejl>/<e<j>. 

<j  >  and,  consequently 

*  <3y>/'<Ey> 

■  <E  >/<E  >. 
a  y 


$ 


In  all  the  calculations  to  be  described  in  Sections  l.l  and  1.2 
the  number  of  cells  was  taken  as  N  =  200,  while  the  averaging  scale 
X  (the  step  of  the  difference  network)  was  assumed  as  h  *  0.1.  A 
perturbation  was  performed  on  the  value  AnQ  in  a  rectangle  consisting 
of  several  cells.  The  shape  of  the  rectangle,  its  position,  and  also 
the  amplitude  and  sign  of  perturbation  were  varied. 

1.1.  It  is  natural  tc  begin  with  a  Coulomb  model  of  a  plasma, 
additionally  assuming  that  T  =  lj  then  v  =  n  ,  a  -  1.  As  calculations 
showed,  the  process  of  the  development  of  ionization  instability  can 
be  arbitrarily  divided  into  two  stages.  In  the  first  stage,  in  a 


certain  characteristic  time  t  *  t*  *  1,  ionization  instability  leads 
to  the  occurrence  of  regular  formations  in  which  n  >  1,  separated 
by  regions  of  lesser  concentration  n  <  .1,  In  shape  and  location 

v 

these  formations  can  be  compared  with  striations  observed  in  the 
experiment.  Figure  2  shows  lines  of  the  level  n®  1  at  time  t  =  1.5 
for  the  case  ft  *  10.  A  perturbation  with  amplitude  0.2  at  t  =  0 
occupied  eight  cells  in  the  central  part  of  the  region  shown  in  Fig.  2. 
The  regions  with  ng  >  1  are  cross-hatched.  The  induced  electric 
currents  are  of  a  vortex  nature  and  flow  along  closed  lines  (Fig.  2). 
Combining  with  the  external  current,  they  cause  the  resulting  current 
to  flow  primari3.y  along  the  striations  within  which  the  intensive 
Joule  heating  is  accompanied  by  ionization.  Attenuations  of  the 
resulting  current  and  Joule  heating,  however,  in  the  "wex Is"'  leads  to 
a  further  decrease  in  them  of  the  concentration  of  free  electrons 
because  of  recombination. 


A  change  In  the  shape  and  position  of  the  perturbation  at  t  *  0 
showed  that  the  shape  and  location  of  the  striations  at  t  'v  t* 
depend  on  the  Initial  conditions.  Thus,  Fig.  3  shows  the  location  of 
the  striations  at  time  t  «  2  when  Ana  *  -0.2.  We  see  that  at  the 
location  of  *the  perturbation  Ang  <  0  a  dip  in  concentration  develops . 
In  general  we  can  note  that  the  positive  sign  of  the  perturbation 
Ang  >  0  Is  caused  by  the  development  of  a  striation  at  the  location 
of  the  perturbation  and,  on  the  other  hand,  the  formation  of  a  "well" 
In  the  case  Ang  <  0.  Let  us  also  note  that  in  all  calculations  the 
characteristic  angle  between  the  average  current  (axis  y)  and  the 
striations  was  approximately  30°-45°. 

In  the  second  stage,  when  t  >  t#  there  is  disruption  of  the 
striations;  these  break  down  into  a  number  of  random  formations 
(•'tongues,"  "loops,"  "bridges,"  etc.).  Figure  4  shows  one  of  the 
moments  of  breakdown  of  the  striations  shown  in  Fig.  2.  The  regions 
with  ng  >  1  are  cross-hatched.  Regular  striations  become,  at  t  >  t#, 
irregular;  in  the  experiments  these  are  compared  with  turbulence. 

A  detailed  description  of  irregular  structure  is  given  in  Section  5. 
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fl  *  10  the  decay  of  the  central  striation  occurs  at  t  2,  while  for 
«  5  it  occurs  at  t  ^  9.5.  If  we  arbitrarily  assume  as  t*  the  moment 
of  decay  of  the  central  striation,  the  dependence  of  t*  on  ft  is  as 
shown  in  Pig.  16.  When  ft  a  ft^  Ci  3  >  the  development  time  t*  tends 
towards  infinity .  This  means  that  when  ft  <  ft  ionization  instability 
does  not  develop  and  no  striation  formation  occurs. 


J2=J  ;  t  si 


Fig.  6. 
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J2  =  $; 


Fig.  10. 
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Let  us  turn  to  an  examination  of  the  average  values.  Figure  17 

shows  the  Hall  ratio  <E  >/<E  >  as  a  function  of  time  t  for  various  ft. 

**  y 

We  should  note  that  the  steady-state  values  B  =  <E  >/<E  >  do  not 

61 I  X  j 

agree  among  themselves  and,  consequently,  do  not  coincide  with  B  . 

V  X 

This  contradicts  the  hypothesis  that  8  -p  >  B.._t  based  on  the  following 

wi<  i  LI 

reasoning.  If  >  8cr,  then  in  the  plasma,  as  a  result  of 

instability,  there  appear  additional  non uniformities  which,  will 

decrease  until  its  value  reaches  8^.  Such  a  foundation  is 

ei  i  cr 

valid,  apparently,  for  small  fluctuations  on  a  uniform  background. 

In  the  developed  stage  of  the  process,  the  fluctuations  are  not  small 

and  the  discussion,  possibly,  becomes  inapplicable.  Figure  18  shows 

the  steady-state  value  <E  >/<E  >  as  t  -*■  «•  as  a  function  of  ft  (solid 

x  y 

line).  From  this  last  graph  we  see  that  there  is  a  tendency  toward 
''saturation*'  of  the  Hall  lelationship  with  an  increase  in  the 
.magnetic  field,  but  the  saturation  itself  occurs  when  ft/fc„0  >  -1-0, 
Finally,  Fig.  19  shows  the  time  variation  of  a  with  a  parametric 
change  in  ft;  Fig.  20  gives  the  dependence  of  the  steady-state  values 
aeff  as  -*■  00  vs.  ft* 
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Pig.  20. 

1.2.  Consideration  of  collisions  of  electrons  with  the  neutral 
component  introduces  certain  qualitative  correction.  Primarily  this 
relates  to  the  behavior  of  n^^  and  t  when  t  >>  t*.  The  value  nmln 
when  t  <*>  now  has  a  certain  limit,  while  <t>,  respectively,  stops 
increasing  without  limit.  This  can  be  seen  in  Pig.  21,  which 
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corresponds  to  the  case  ft  =  10  when  =  0.2.  The  ratio  of  the 
initial  temperature  Tg  to  the  ionization  potential  I  is  taken  as 
kn  *  0.1.  The  ratio  Aj/A^  *  1,  A^  -  0.  The  lines  of  level  ng  *  1  and 
the  field  of  electric  current  densities  at  moment  of  time  t  *  1,  5* 
and  20  are  shown  in  Pigs.  22-2H. 

Calculations  were  also  conducted  for  other  values  of  ft.  It  is 
characteristic  that  when  taking  into  account  collisions  with 
neutrals,  the  effect  of  "saturation”  of  the  Hall  ratio  with  increasing 
ft  becomes  more  clearly  expressed  (Fig.  18,  dashed  line). 


t 


j}*/0 ;  tm3 
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For  an  explanation  of  the  influence  of  degree  of  ionization  of 
the  plasma  on  the  nature  of  the  development  of  ionization  instability 
we  perform  calculations  for  various  ratios  of  Ag/A-^.  Figures  25-28 
'show  the  results  of  a  numerical  solution  of  a  version  analogous  to 
that  examined  above  but  with  Ag/A^  3  10,  i.e.,  in  the  case  of  a 
weakly-ionized  plasma.  We  can  note  that  the  process  of  development 
of  strlations  occurs  somewhat  more  slowly  than  when  Aj/'A^  =  1.  The 
qualitative  dependence  of  t*  on  Ag/A^  w^en  ft  *  10  is  shown  in 
Pig.  29. 
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2.  Let  us  change  the  statement  of  the  problem  as  follows.  Let 
us  assume  that  the  dielectric  walls  are  spaced  farther  apart,  while 
the  ideally  sectioned  electrodes  ab  **  1  and  cd  *  1  are  located  right 
in  the  middle.  Sections  Aa,  bB,  Cc,  and  dD  are  ideal  dielectrics. 

The  current  density  on  the  electrodes  of  the  process  is  given  as 
j  *  -1.  The  number  of  cells  N  *  600  (h  “  0.1).  Otherwise',  the 

J 

statement  is  as  in  section  1.  Let  us  note  that  there  is  as  yet  no 
need  to  Introduce  perturbation  into  the  initial  conditions,  since 
there  is  no  trivial  solution. 

2.1.  Coulomb  plasma.  Figures  30-3^  show  the  lines  of  the  level 
and  the  vector  field  of  the  electric  current  density  for  ft  *  10  and 
times  t  *  1,  1.8,  3,  7*  and  18.2. 

At  the  Initial  moment  of  time  the  electron  concentration  is 
uniform;  therefore  the  electric  current  field  is  perturbed  only  by 
the  Hall  effect.  Near  the  electrodes  the  current  density  is  highest; 
here  regions  of  increased  concentration  ne  are  formed.  Consequently, 
the  effective  Hall  parameter  here  abruptly  decreases,  while  the 
currents  "approach"  In  direction  the  vector  of  the  average  electric 
field  and  flow  primarily  along  the  electrodes.  Powerful  current 
layers  are  formed  which  depart  from  the  electrodes  and  propagate  along 
the  dielectric  portions  Aa  and  Cc  (Fig.  30).  Within  the  layers  there 
is  Intense  plasma  ionization.  Near  the  points  b  and  c  the  electron 
concentration  rapidly  increases  and  reaches  a  limiting  value. 

H2 
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Because  a  Is  finite,  the  currents  from  the  layers  are  diffused  and 
directed  toward  opposite  electrodes. 
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Fig.  31*. 

After  a  certain  time  (Fig.  31),  between  the  electru^es  there 
appears  a  region  of  increased  ionization  (striation)  along  which  the 
electric  current  contracts  and  there  is  intense  release  of  Joule 
heat,  causirfg  an  increase  in  concentration  ne.  The  maximum  concentra¬ 
tions  ne  about  points  b  and  c  decreases  somewhat,  since  an  ever 
increasing  considerable  part  of  the  curr'ent  tends  from  the  electrodes 
along  the  striation.  This  drop  in  n_Q_  is  shown  in  Fig.  35,  where, 

nicLjv 

in  addition  to  nmax  and  nmin»  there  are  depicted  the  average  and 
effective  values  as  functions  of  time.  The  point  of  maximum 
concentration  appears  to  lie  within  the  striation. 

In  the  second  stage,  when  t  >  t*  (for  the  given  variant  t*  5), 
the  striation  breaks  down  into  several  current  pinches  located  almost 
vertical  from  electrode  to  electrode.  Approximately  the  same  picture 
is  observed  during  experimental  study  of  striations  on  finitely 
sectioned  electrodes. 

As  before,  the  transition  from  the  second  stage  is  accompanied  by 
the  establishment  of  average  values  (excluding  <t>  md  nmin),  despite 
the  absence  of  stable  stationary  distribution  of  n  and  the  electric 
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current  field.  The  asymptotic  value  of  <Ex>/<Ey>  at  t  +  “  became  equal 
to  3.3  for  fi  *  10  Instead  of  the  value  of  5.5  given  in  Section  2,  while 
°eff  ^his  case  tencte  towards  0.35. 


Calculations  performed  for  various  values  of  ft  showed  that  the 
appearance  of  ionization  instability  begins  at  ftcr  y  3. 

2.2.  Consideration  of  collisions  of  electrons  with  neutrals, 
just  as  in  the  problem  of  Section  1,2,  leads  to  asymptotic  establish¬ 
ment  of  <t>  and  njnin .  We  can  also  again  point  out  the  slower 
transition  to  the  second  stage  of  ''pinching"  of  the  currents  and 
decay  of  the  strlations.  The  results  of  numerical  solution  of  the 
problem  of  Section  2.1,  with  consideration  of  collisions  of  elt  .  ?ns 
and  neutrals  at  n  =  10  and  Ag/A1  =  1,  are  given  in  Pigs.  36-38  for 
times  t  *  1,  5,  and  8. 
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The  time  variation  of  effective  and  average  values  is  shown 
in  Fig.  39.  Let  us  note  that  the  dip  in  n__„  occurring  during  the 

luaX 

striation  formation  process  then  disappears,  unlike  the  case  of 
,  Coulomb  plasma. 

1  3*  Let  us  examine  the  interaction  of  nonequilibrium  ionized 

plasma  with  solid  ideally  conducting  electrodes.  Here  we  limited 
ourselves  to  the  case  of  Coulomb  plasma  a  -  1. 

3.1.  First  let  us  examine  a  statement  of  the  problem  similar  to 

that  in  Section  1.  At  the  initial  moment  t  =  0,  j  ,  j  ,  and  n  are 

y  “ 

uniformly  distributed  in  space  with  values  i  *  0,  J  *  -1,  and 

x  y 

n  *  1.  An  external  magnetic  field  is  applied  to  the  entire  volume 
1  ABCD  (AB  *  CD  *  1,  AD  =  BC  *  2),  and  its  value  is  such  that  at  the 
'initial  moment  cf  time  the  Hall  parameter  ftQ  =  10.  The  horizontal 
section  AB  and  CD  of  the  boundary  of  the  region  (see  Fig.  1)  are  solid 
ideally  conducting  electrodes  instead  of  ideally  sectioned  electrodes 
examined  in  Section  1. 

Figures  40-43  show  the  three-dimensional  distribution  of  the 
electric  current  density  components  Jx,  jy  and  electron  concentration 
n  ,  respectively,  for  times  t  *  0.2,  0.4,  1.14,  and  4.0. 

C 

Let  us  remember  that  for  the  problem  examined  in  Section  1  the 

initial  data  j  =  0,  j  *  -1,  n  «  1  are  the  steady-state  solution, 

j  x  y  “ 

j  In  the  problem  with  solid  electrodes  such  values  of  j  .  j  ,  and  n^ 

I  x  y  e 

1  satisfy  the  equations,  but  do  .not  satisfy  the  boundary  conditions  on 

1  the  solid  electrodes:  J  =  -ft.i  ,  Therefore  here  there  is  no  need 

•  x  y 

i  to  introduce  perturbations  into  the  initial  data  in  order  to 

Initiate  the  nonstationary  process.  In  the  initial  stage,  because 
of  the  strong  Hall  effect  (ft  =  10),  the  current  is  concentrated  on  the 
left  end  of  the  upper  electrode  AB  (point  A  =  a)  and  on  the  right  end 
of  the  lower  electrode  CD  (at  point  C  *  c).  Corresponding  to  this, 
in  the  upper  left  and  lower  right  corners  of  the  region  there  is 
localization  of  Joule  heating,  which  leads  to  increased  ionization  cf 
plasma  and  the  achievement,  at  t  =  0.2,  of  maximum  electron  concentration 

I 
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The  time  variation  of  effective  and  average  values  is  shown 
in  Pig.  39.  Let  us  note  that  the  dip  in  nmgx  occurring  during  the 
striation  formation  process  then  disappears,  unlike  the  case  of 
Coulomb  plasma. 

3.  Let  us  examine  the  interaction  of  nonequi librium  ionized 
plasma  with  solid  ideally  conducting  electrodes.  Here  we  limited 
ourselves  to  the  case  of  Coulomb  plasma  a  =  1. 


3.1.  First  let  us  examine  a  statement  of  the  problem  similar  to 

that  in  Section  1.  At  the  initial  moment  t  =  0,  j  ,  j  ,  and  n  are 

^  y  s 

uniformly  distributed  in  space  with  values  i  »  0,  j  =  -1,  and 

x  y 

n  *  1.  An  external  magnetic  field  is  applied  to  the  entire  volume 
ABCD  (AB  *  CD  »  1,  AD  *  BC  *  2),  and  its  value  is  such  that  at  the 
initial  moment  cf  time  the  Hall  parameter  ftQ  =  10.  The  horizontal 
section  AB  and  CD  of  the  boundary  of  the  region  (see  Pig.  1)  are  solid 
ideally  conducting  electrodes  instead  of  ideally  sectioned  electrodes 
examined  in  Section  1. 


Figures  40-43  show  the  three-dimensional  distribution  of  the 
electric  current  density  components  Jx,  and  electron  concentration 
n&,  respectively,  for  times  t  «  0.2,  0.4,  1,14,  and  4.0. 


Let  us  remember  that  for  the  problem  examined  in  Section  1  the 
initial  data  *  0,  j  *  -1,  nft  *  1  are  the  steady-state  solution. 
In  the  problem  with  solid  electrodes  such  values  of  j  ,  j 
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satisfy  the  equations,  but  do  ..not  satisfy  the  boundary  conditions  on 

the  solid  electrodes:  j  =  -ft,*  .  Therefore  here  there  Is  no  need 

x  y 

to  Introduce  perturbations  into  the  initial  data  In  order  to 
initiate  the  nonstationary  process.  In  the  initial  stage,  because 
of  the  strong  Hall  effect  (ft  =  10),  the  current  is  concentrated  on  the 
left  end  of  the  upper  electrode  AB  (point  A  =  a)  and  on  the  right  end 
of  the  lower  electrode  CD  (at  point  C  »  c).  Corresponding  to  this. 

In  the  upper  left  and  lower  right  corners  of  the  region  there  is 
localization  of  Joule  heating,  which  leads  to  increased  ionization  cf 
plasma  and  the  achievement,  at  t  =  0.2,  of  maximum  electron  concentration 
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ng  a  3  (Pig.  40).  Beginning  at  moment  t  %  1  (Pig.  4l)  the  regions  of 

increased  electron  concentration  (ng  >  1)  are  propagated  along  the 

vertical  sections  of  boundaries  AD  and  CB,  forming  along  them  layers 

of  increased  concentration  n  ,  while  in  the  center  of  the  region  the 

electron  density  drops.  The  current  passes  primarily  along  these 

layers,  passing  between  them  and  agitating  the  plasma  in  the  center 

of  the  region  (Pig.  42).  By  time  t  =  4  (Fig.  43)  the  layer  of 

increased  electron  concentration  has  extended  along  the  entire 

boundary  of  the  region,  while  inside  on  the  background  of  electron 

concentration  ne  %  0.5>  the  current  density  is  low  in  comparison  with 

the  density  on  the  boundary,  and  is  randomly  distributed  in  space. 

Thus,  the  local  nature  of  distribution  j  ,  j  ,  and  n  in  space  differs 

a  y  v 

substantially  from  the  nature  of  the  distribution  with  ideally 

sectioned  electrodes  (see-  Figs.  2-4).  This  leads  to  a  substantial 

difference  in  the  integral  characteristics  (cf.  Pig.  5  and  Pig.  44). 

The  establishment  of  magnitudes  <Ex>/<Ey>,  and  with  time  here 

occurs  by  the  time  t  *  1,  l.e.,  considerably  sooner  than  in  the  case 

of  ideally  sectioned  electrodes.  Turbulent  resistance  is  greater  by 

a  factor  of  2.5,  while  the  ratio  <E  >/<E  >  is  an  order  of  magnitude 

x  y 

smaller. 
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n  -  3  (Pig.  40).  Beginning  at  moment  t  %  1  (Pig.  41)  the  regions  of 

increased  electron  concentration  (ng  >  1)  are  propagated  along  the 

vertical  sections  of  boundaries  AD  and  CB,  forming  along  them  layers 

of  increased  concentration  ne,  while  in  the  center  of  the  region  the 

electron  density  drops.  The  current  passes  primarily  along  these 

layers,  passing  between  them  and  agitating  the  plasma  in  the  center 

of  the  region  (Pig.  42).  By  time  t  =  4  (Fig.  43)  the  layer  of 

increased  electron  concentration  has  extended  along  the  entire 

boundary  of  the  region,  while  inside  on  the  background  of  electron 

concentration  ng  ^  0.5,  the  current  density  is  low  in  comparison  with 

the  density  on  the  boundary,  and  is  randomly  distributed  in  space. 

Thus,  the  local  nature  of  distribution  j  ,  j  ,  and  n  in  space  differs 

a  y  “ 

substantially  from  the  nature  of  the  distribution  with  ideally 
sectioned  electrodes  (see-  Pigs.  2-4).  This  leads  to  a  substantial 
difference  in  the  integral  characteristics  (cf.  Pig.  5  and  Pig.  44). 
The  establishment  of  magnitudes  <Ex>/<Ey>,  and  creflf.  with  time  here 
occurs  by  the  time  t  *  1,  l.e.,  considerably  sooner  than  in  the  case 
of  ideally  sectioned  electrodes.  Turbulent  resistance  is  greater  by 
a  factor  of  2»5,  while  the  ratio  «E  >/<E  >  is  an  order  of  magnitude 

A  J 

smaller. 
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3.2.  Now  let  us  examine  the  same  problem  but  in  more  complex 
geometry.  We  will  consider  that  an  electric  current  flows  between 
the  solid  ideal  electrodes;  these  comprise  only  a  part  of  ab  and  cd 
of  the  horizontal  sections  of  boundaries  ab  and  cd  (see  Pig.  1). 

An  analogous  problem  was  examined  in  Section  2,  but  with  ideally 
sectioned  electrodes;  therefore  henceforth  the  results  will  be 
compared  with  those  of  Section  2.  By  time  t  *  1  (cf.  Pig.  30  and 
Pig.  45)  which  is  characteristic  for  both  problems  we  have  a  boundary 
layer  of  increased  concentration  of  electrodes,  which  was  discussed 
in  Section  2.  Let  us  note  that  because  of  the  stronger  concentration 
of  current  at  the  ends  of  the  solid  electrodes  a  and  c  and  on  the 
ends  of  the  ideally  sectioned  electrodes,  the  formation  of  this 
layer  occurs  more  intensely.  In  the  center  of  the  region  at  this 
moment  with  solid  electrodes  the  current  has  a  primarily  horizontal 
direction,  while  with  the  ideally  sectioned  electrodes  it  was 
vertical.  The  three-dimensional  distribution  of  current  density 
and  j  and  electron  concentration  ng  are  shown  in  Pigs.  45-49, 
respectively,  for  times  t  *  1,  3,  5,  7,  and  10. 


At  time  t  *  3  the  current  in  the  center  of  the  region,  having  a 
horizontal  direction  (Fig.  46),  forms  a  striation  with  increased 
electron  concentration.  By  time  t  *  5  the  striation  expands  and 
disintegrates  into  individual  horizontal  "streaks.”  / 

/ 

This  is  explained  by  the  fact  that  E  <<  E  between  the  electrodes 
because  of  the  influence  of  the  boundary  condition  Ex  »  0  on  the 
electrodes  themselves.  Since  the  actual  Hall  parameter  is  large  here 
(since  concentration  ng  is,  on  the  average,  sm'ail),  the  electric 
current  flows  almost  perpendicular  to  the  electric  field  which  can 
effectively  be  identified  with  Eg.  Thus,  the  electric  current  in  the 
interelectrode  region  flows  primarily  parallel  to  the  ideal  electrodes. 
Figures  48  (t  *  7)  and  4,9  (t  «  10)  show  further  agitation  of  the 
plasma;  a  characteristic  feature  is  that  the  current  breaks  down  into 
.individual  streaks,  having  primarily  horizontal  direction. 

Let  us  note  that  the  flow  of  current  in  ideally  sectioned 
electrodes  with  the  same  geometry  and  the  same  parameters  ha*  a 
different  nature  (see  Figs.  30-34).  There,  by  time  t  *  1.8,  there  is 
fonmed  a  bundle  with  increased  electron  concentration  (striation) 
directed  from  the  right  end  of  the  upper  electrode  (point  b)  to  the 
left  end  of  the  lower  electrode  (point  c).  Then  the  plasma  is 
agitated,  but  the  electric  current  remains  concentrated  in  the 
interelectrode  space  and  breaks  down  into  streaks  which  have 
primarily  vertical  direction  and  which  are  connected  by  bridges 

Such  a  difference  in  the  passage  of  the  electric  current  leads 

to  varying  changes  in  the  average  values  with  time  (cf.  Fig.  35  and 

?ig.  50).  Turbulent  resistance  of  the  plasma  is  approximately  twice 

as  great,  while  the  ratio  <E  >/<£  >  2l  0.5  instead  of  the  ratio  ~*2  in 

a  y 

the  case  of  ideally  sectioned  electrodes. 

4.  Let  us  examine  the  Influence  of  nonuni  fcrmity  of  me 
magnetic  field  on  the  passage  of  current  between  ideally  sectioned 
electrodes  in  the  geometry  of  Section  2.  Let  us  assume  that  the 
geometry  of  the  problem  and  the  initial  conditions  at  initial  moment 
of  time  t  *  0  are  analogous  to  those  In  Section  2.  The  problem  differs 
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from  that  in  Section  2  in  that  the  external  magnetic  field  occurs  only 

between  electrode* j  in  the  remaining  sections  of  the  regions  it  is 

absent-  Along  lines  ad  and  be  there  is  an  abrupt  drop  in  the 

magnetic  field  from  H  «  10  to  H  *  0.  The  distribution  of  j  ,  j  ,  and 

x  y 

n  is  shown  in  Pigs*  51 >  52,  and  53,  for  times  t  ®  1,  5?2,  and  11.6, 

V 

respectively.  In  comparison  with  the  corresponding  problem  of 
Section  2  (see  Figs.  30-31*)  the  picture  of  j  ,  and  ne  differs 
substantially.  There  are  no  layers  of  increased  concentration  along 
the  dielectric  wall  (since  ft  in  AadD  and  bBCc).  In  these  regions 
there  is  "barrel-shaped"  flow  of  current  from  the  upper  electrode  to 
the  lower  one.  Discontinuity  of  the  magnetic  field  along  ad  and  be 
causes  discontinuity  in  ft,  which  leads  to  discontinuity  of  the 
electric  field  intensity  and  the  electric"  current  density  along  ad 
and  be.  With  time,  because  of  instability,  there  appears  turbulent 
resistance  of  the  plasma  in  the  space  between  the  electrodes; 
therefore  the  electric  current  flows  from  electrode  to  electrode 
primarily  along  regions  AadD  and  bBCc  (Figs.  52-53)* 


Fig.  51. 
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The  maximum  currents  are  concentrated  near  the  magnetic  field 
and  electrode  discontinuities,  causing  maximum  ionization  there;  in 
the  remaining  space  there  is  no  current.  Figure  5^  shows  the  integral 
values.  A  comparison  of  them  with  integral  values  of  the  corresponding 
problem  in  Section  2  (Fig.  35)  shows  that  the  end  effects,  due  to  the 
cutoff  of  the  magnetic  field,  leads  to  a  drop  (two-fold)  in  the  Hall 
intensity  of  the  electric  field  and  an  increase  in  turbulent 
resistance. 

5.  The  problem  of  ionization  instability  in  an  unbounded  volume, 
as  can  be  seen  from  the  initial  equations,  has  no  characteristic 
dimension.  However,  when  solving  it  numerically,  the  characteristic 
dimension  is  the  scale  of  averaging  equal  to  the  step  of  the  difference 
net  X  =  h.  When  solving  the  problem,  taking  into  account  diffusion 
processes  (radiation,  thermal  conductivity,  etc.),  the  characteristic 
dimension  of  the  problem  will  be  the  diffusion  length  of  these 
processes.  From  this  point  of  view,  the  step  of  the  difference  net 
can  be  treated  as  the  diffusion  length  of  certain  physical  processes. 


while  a  change  in  step  of  the  net  corresponds  to  a  change  In  certain 

physical  parameters  (mean- free  path  of  radiation,  coefficient  of 

thermal  conductivity,  etc.}.  Let  us  trace  the  influence  of  the  scale 

of  averaging  on  the  development  of  ionisation  instability.  In  Section 

1  we  examined  the  problem  for  a  number  of  computation  cells  N  *  200 

(h  *  h  *  0.1),  and  a  perturbation  was  introduced  into  eight  cells  in 
a  y 

the  center  of  the  region.  Let  us  examine  this  same  problem  for 

N  *  800  (h_  *  h„  *  0.05)  with  a  perturbation  in  32  cells.  A  decrease 
a  y 

in  the  step  of  the  difference  net  leads,  at  time  t  «  1.5,  to  a  somewhat 
different  location  of  the  striation  (Pig.  55)  than  in  Pig.  2.  At  the 
location  of  the  same  disturbance  there  occur  two  striations  whose 
thicknesses,  as  before,  are  equal  to  several  scales  of  averaging 
(several'  “steps  of  the  difference  net).  Near  the  boundary  the 
striation  picture  is  disrupted  slightly,  which  is  explained  by  the 
influence  of  boundary  conditions.  The  formation  of  striations  in  this 
case  occurs  more  intensely  and,  by  moment  t  *  1.5,  we  find  indications 
of  the  appearance  of  an  irregular  picture  which,  by  moment  t  *  7, 
has  little  in  common  with  the  picture  in  Pig.  4.  The  development  of 
instability  at  N  *  800  from  a  perturbation  occupying  four  cells  in 
the  center  of  the  region  is  shown  Pigs.  56-58,  corresponding  to 
moments  of  tiime  t  «  1.4,  4.6,  and  19.0.  The  local  picture  of  Jx,  jy, 
and  n  differs  from  the  previous  problem.  By  time  t  *  1.4  (Fig.  56) 
there  are  five  striations;  the  central  striation  is  located  at  an 
angle  close  to  ir/3  to  the  total  current.  Then  the  striations  are 
disrupted  (Pig.  57)  and  by  moment  t  **  19  (Pig.  58)  there  is  strongly 
developed  Instability  in  the  plasma.  Let  us  note  that  with  passage 
along  such  a  plasma  the  current  breaks  down  into  elongated  vertical 
streaks  (5-10  characteristic  dimensions)  which  are  connected  by 
short  bridges  (1-2  characteristic  dimensions) .  Let  us  examine  the 
change  in  time  of  the  average  values  which,  for  N  »  8C0,  are  shown 
in  Pig.  59.  Despite  the  varying  distribution  of  electron  concentra¬ 
tion  with  time  during  complete  perturbation  in  the  32nd  and  8th  cells 
the  integral  characteristics,  beginning  with  moment  t*  *  2  coincide 
with  good  accuracy. 

With  a  narrower  step  the  drop  in  ceff  and  0eff  =  <Ex>/<Ey> 
occurs  more  intensely;  thus,  oeff  decreases  from  6  at  N  *  200  to  3*2 
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at  N  «  800.  This  is  evidently  associated  with  the  drop  in  $  .  with 

decreasing  step  of  the  net.  efi 


Q-&  l*li 


Pig.  55. 


Pig.  56 


Let  us  illustrate  the  influence  of  a  change  in  characteristic 
length  of  averaging  X  in  yet  another  problem.  Let  us  assume  that  the 
ideally  sectioned  electrodes  occupies  the  central  region  (10  intervals, 
h  *  0.1)  of  the  horizontal  sections  of  the  boundary  along  which  there 
are  2d  calculation  intervals.  The  vertical  walls  occupy  20  intervals. 

The  geometry  of  this  problem  is  the  same  as  that  of  Section  2. 
Again  we  will  consider  the  plasma  to  be  Coulomb  and,  when  t  »  0, 
j  *  0,  j  *  -1  in  the  space  between  the  electrodes „  while  ne  *  1 
throughout  the  region.  The  distribution  of  Jx,  jy,  and  ne  at  moments 
t  ■  3  and  5  is  shown  in  Figs.  60  and  6l  (only  the  left  halves  of  the 
figures  are  shown).  The  passage  of  current  in  this  problem  repeats 
all  the  features  noted  in  Section  2  (see  Figs.  30~3*O  • 

The  distribution  of  jx<,  jy,  and  n  in  this  same  geometry,  but  with 
h  *  -0.05,  at  moments  t  *  3,  5,  and  7  are  shown  in  Figs.  62-64.  A 
comparison  of  Figs.  60  and  62,  and  also  Si  and  63,  snows  that 
ionization  instability  develops  more  intensely  with  a  shorter  length 
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of  averaging,  while  the  turbulent  stage  is  achieved  by  the  plasma 
earlier.  There  is  formations  of  current  streaks  having  width 
'''A  *  h,  elongated  in  the  vertical  direction  (M.GX),  connected  by  short 
('vX  -  2X)  bridges.  This  can  be  seen  In  Pig.  6*1,  which  shows  the 
distribution  of  Jv,  and  n  when  t  =  7. 


60. 


61. 


3 
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Pig.  6H. 


The  conclusion  as  to  the  growth  of  turbulent  resistance  oeff 

and  a  drop  in  ratio  <E  >/<E  >  with  decreasing  X  is  confirmed  by 

y 

Pig.  65,  which  shows  the  change  in  integral  values  with  time 
(X*h»0.05-  solid  lines,  X  *  h  *  0.1  —  dashed  lines). 
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i  'v  <n>X  the  current  should  "leave4  the  streak,  "breaking”  it.  Such 
"broken"  streaks  can  cover  the  region  only  by  forming  "large-scale 
layers"  (cf.  [18])  of  thickness  'vJt.  Along  the  continuous  boundary 
of  the  large  layers,  electric  currents  pass  in  turn.  As  numerical 
calculations  show,  here  there  are  formed  current  bundles  with  increased 
concentration  ne,  which  can  be  identified  with  striations.  The 
striations  have  a  thickness  of  the  order  of  several  scales  of  averaging 
and  are  elongated  (parallel  to  each  other)  at  a  certain  angle  to  the 
average  current  (see  Pig.  58).  The  value  of  the  angle  corresponds, 
apparently,  with  the  direction  of  most  intense  release  of  Joule  heat. 
This  corresponds  with  the  fact  that  electron  concentration  and 
current  density  in  the  striations  are  considerably  higher  than  in 
the  streaks.  In  the  calculations  we  observed  the  angle  -  j.  The 
value  ir/4  corresponds  to  the  solution  of  the  problem  on  the  develop¬ 
ment  of  one-dimensional  perturbations. 

7.  Let  us  present  the  results  of  characteristic  calculations 
with  H  »  1,  described  above,  in  the  form  of  Table  1  (Pig.  66).  Let 
us  explain  the  designations.  The  region  occupied  by  the  plasma  Is 
arbitrarily  divided  into  two  types:  The  first  type  Includes  a  region 
whose  horizontal  sections  are  entirely  covered  by  the  electrodes; 
the  second  type  Includes  the  region  whose  horizontal  sections  contain 
dielectric  portions.  The  other  designations  are  clear  from  what  has 
been  said  above.  The  values  of  the  effective  parameters,  presented 
in  the  table,  satisfy  the  following  relationship  quite  accurately: 

The  values  of  che  effective  parameters,  presented  in  the  table, 
satisfy  the  following  relationship  quite  accurately: 


°eff 


eff 


'OIT 


(5.1) 


where  fln  and  cn  are  values  corresponding  to  uniform  distribution  of 

U  U 

the  plasma  parameters,  i.e. ,  Cq  *  1,  9.  m  —  «  H.  Let  us  turn  our 
attention  to  the  fact  that  relationship  (5.1)  is  invariant  with  respect 
to  type  of  electrode  and  region,  scale  of  averaging  X,  and  also  ratio 
A2^A19  which  characterizes  the  presence  of  a  neutral  component.  The 
validity  of  (5.1)  is  also  confirmed  by  the  results  of  all  other 
numerical  calculations.  An  exception  Is  that  line  of  the  table. 
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Indicated  by  an  asterisk,  which  coi’responds  to  the  case  when  the 
magnetic  field  occurs  only  in  interelectrode  space.  Let  us  remember,  j 
however,  that  the  electric  current  in  this  case  flows  almost  completely  j 
from  the  .Interelectrode  region  and  occurs  laminarly  along  a  non-  | 

magnetized  plasma,  while  (5«1)  obtains  for  a  plasma  with  a  strongly  j 

developed  Instability.  ' 


Table  1. 


Points 

Dimension 

of 

averaging 

Type 

of 

region 

Type 

_  £t 

Of 

electrode 

eeff 

aeff 

ft0 

<n> 

<g> 

Remarks 

1. 

0.1 

1 

ideally 

sectioned 

3.5 

B 

10 

B 

11.2 

a  neutral  com¬ 
ponent  exists 

2. 

0.1 

1 

ideally 

sectioned 

5.5 

r 

0.55 

10 

g 

17 

3. 

0.1 

1 

solid 

2 

0.2 

10 

B 

14.8 

4. 

0.05 

1 

ideally 

sectioned 

2.5 

0.32 

10 

B 

14.4 

5: 

0.1 

2 

solid 

ESI 

0.06 

10 

B 

B 

6. 

0.05 

2 

ideally 

sectioned 

2 

0.2 

10 

B 

B 

1 

0.1 

1 

ideally 

sectioned 

3.6 

0.75 

5 

B 

B 

1 

•8. 

0.1 

2 

ideally 

sectioned 

fl 

0.3 

10 

— 

— 

magnetic  field 
between  the 
electrodes 

Relationship  (5*1)  can  be  obtained  approximately  as  follows. 

Analogous  to  [18],  let  us  examine  a  layer  of  plasma  between  strlations, 

developed  along  the  streaks.  Let  us  direct  axis  x  perpendicular  to 

the  streaks  and  axis  y  along  them  (parallel  to  the  average  current). 

Setting  <8>  »  1  we  can  consider  that 

streaks,  J  »  j, 

j  . 

the  equation 


* 

O  i 


xy  n 


x  -\fy‘ 

(obviously),  while  between  them,  Q  >>  1. 


Actually ,  in  the 


Thus  in 


J*  +  njy  ' 

the  component  j  can  be  discarded,  i.e.. 
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4  ,  2e  *  neEx 

“x  ftEx  H  * 

Let  us  average  this  equality,  without  considering  the  correlation 
between  n  and  E  : 

t-  A 

'  <V<V 

<i  >  s  — ...  ^  _ 2 _ £_ 

Jy  H  -  H 


From  this, 


<j-  >  <n  >  <E  >  n 
t  «  -X  »  e  -X  =  :_e 

eff  <E  >  H  <E  >  H  peff  * 

J  J 


Since  in  numerical  calculations  <ne>  *  1,  we  can  write,  finally, 
setting,  nQ  «  aQ  «  1, 


aeff  “  °0  aQ  * 


In  the  derivations  we  did  not  consider  currents  in  striations 
on  the  boundaries  of  the  layers.  But  since  the  thickness  of  each 
stratum  is  smaller  than  the  thickness  of  the  layer  by  a  factor  of 
while  the  current  in  it  is  limited,  then  when  »  1  the 
w contribution"  of  striations  to  the  average  parameters  of  the  plasma 
can  be  disregarded.  Finally,  if  we  assume  that  when  A  <<  1  the 
equality  *  Bcr  is  nonetheless  valid,  we  get 


°eff  *  «0Pcr* 


We  should,  however,  stress  that  such  an  approach  must  be  founded. 
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CHAPTER  6 
CONCLUSIONS 

On  the  basis  of  an  analysis  of  the  results  of  numerical  solution 
of  the  problem  of  the  development  of  ionization  Instability  in  a 
magnetized  plasma  we  can  draw  the  following  conclusions. 

1.  The  physical  model  used  qualitatively  reflects  the  essence 
of  the  phenomenon  and  the  characteristic  moment  known  from  the 
experiment : 

a.  The  nonlinear  mechanism  of  Joule  heating  and  transfer  of 
energy  during  collisions  of  electrons  with  a  heavy  component  causes, 
in  characteristic  time  t»,  the  appearance  of  ionization  instability 
as  regular  formations  —  striations. 

b.  In  later  stages  cf  the  process  (t  >>  t*)  there  is  decay  of 

the  primary  striations,  leading  to  irregular  three-dimensional 
distribution  of  electron  concentration  and  electric  current  density. 
Nonetheless,  when  t  »  t*  we  speak  of  the  structure  of  an  Irregular 
picture,  which  consists  of  the  following.  An  electric  current  passes 
mainly  in  the  "streaks'*  parallel  to  the  direction  of  the  average 
current.  The  thickness  of  the  streaks  and  the  difference  between 
them  'Wi,  while  their  length  where  X  is  the  characteristic 

scale  of  averaging.  Between  the  streaks  the  current  density  and 
concentration  nfi  are  relatively  low.  The  streaks  are  joined  into 
layers  along  whose  continuous  boundaries  there  are  strong  electrical 
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currents,  forming  striations  with  high  electron  concentration.  The 
angle  between  striation  and  the  direction  of  the  average  current 

TT  If 

^  the  thickness  of  the  striations  is  of  the  order  of  several 
scales  of  averaging. 

c.  On  the  electrodes  and  the  dielectric  regions  adjacent  to 
them  there  are  layers  with  increased  electron  concentration;  these 
were  observed  by  Kerrebrock. 

d.  The  presence  of  a  nonmagnetized  region  of  plasma  leads  to 
emergence  of  the  electrical  field  into  this  region. 

2.  Despite  the  locally  nonstationary  picture  of  developed 
instability  there  occurs  an  establishment  of  average  and  effective 
parameters  of  the  plasma  in  certain  characteristic  times  t*.  Here  we 
should  note  the  following: 

a.  .  The  steady-state  average  and  effective  parameters  of  the 
plasma „  generally  speaking,  depend  on  the  type  of  electrode,  the 
geometry  of  the  region,  the  presence  of  a  neutral  component,  and  the 
scale  of  averaging  A,  and  do  not  depend  on  the  initial  condition. 

b.  The  effective  parameters  of  the  plasma  are  connected  by  the 
relationship 

aeff  *  fl^eff 

which  is  invariant  with  respect  to  the  geometry  of  the  region,  the 
type  of  electrode,  the  presence  of  a  neutral  component,  and  the 
"Seale  of  averaging  A. 

3.  We  have  refined  the  statistical  picture  of  strongly  developed 
instability: 

a.  There  is  angular  anisotropy  in  the  distribution  of  the 
electron  concentration  relative  to  the  direction  of  the  average 
current,  which  is  released. 


b.  There  is,  generally  speaking,  no  mirror  symmetry  of  distribu¬ 
tion  of  electron  concentration  with  area. 

b .  The  boundary  conditions  have  a  substantial  Influence  on  the 
process  of  development  of  instability  and  the  final  state  of  the 
plasma: 

a.  Between  the  ideally-sectioned  electrodes  the  current  passes 
primarily  in  a  direction  orthogonal  to  the  plane  of  electrodes. 

b.  Between  the  solid  ideal  electrodes  the  current  in  the 
interelectrode  region  is  directed  primarily  parallel  to  the  electrode 
plane. 


c.  The  average  and  effective  parameters  of  the  plasma,  with 
t  »  t*,  as  already  noted,  depend  on  the  boundary  conditions. 

5.  A  decrease  in  the  characteristic  scale  of  averaging  leads  to 

a  reduction  in  the  first  stage  of  the  process  (time  t*  of  the 

formation  of  striations  and  decrease  of  3eff  (however  3eff  is  always 

greater  than  3„„)). 

cr 

6.  Consideration  of  collisions  of  electrons  with  a  neutral 
component  does  not  change  substantially  the  qualitative  nature  of  . 
the  development  of  instability  and  the  final  state  of  the  plasma. 

The  authors  would  like  to  thank  A.  A.  Samarskiy  for  his  valuable 
discussions  and  Ye.  Ye.  Myshetskaya  for  compiling  t/ne  programs 
and  carrying  out  the  calculations. 
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